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I. INTRODUCTION 
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Isospin singlet (pn) pairing as well as quartetting in nuclei is expected to arise near the symmetry 
line N = Z. Empirical values can be deduced from the nuclear binding energies applying special 
filters. Within the local density approximation, theoretical estimates for finite nuclei are obtained 
from results for the condensation energy of asymmetric nuclear matter. It is shown that the isospin 
singlet condensation energy drops down abruptly for \N — Z\ ~ 4 for medium nuclei in the region 
ON ' A — 40. Furthermore, a-like quartetting and the influence of excitations are discussed. 

Keywords: proton-neutron pairing, superfluidity, strongly coupled systems, nuclear matter, a- 
<**] ■ particle matter. 
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~^J ■ The structure of nuclei is understood to a large extent within a single-particle approach, the nuclear shell model 

which was introduced in nuclear theory 50 years ago M . Its microscopic foundation should be given from a quantum 

statistical treatment of the many-nucleon system, where in general also the concept of temperature can be introduced 

to describe excited nuclei in a dense medium. Within the mean-field approximation, the single-particle approach can 

ON . be obtained . 

However, due to the nucleon-nucleon interaction also correlations occur which cannot be accounted for within this 
quasi-particle approach. Two-particle correlations and higher order correlations as precursors of corresponding bound 
states may become of importance in regions where the nucleon density is low. A systematic treatment of the effect of 
correlations can be given in terms of the spectral function, see, e.g., pM and further references given therein. 

Free two-nucleon states are described by scattering phase shifts or, more generally, by the T matrix. In the isospin 
singlet (S = 1, T = 0) channel where the interaction is stronger than in the isospin triplet (S = 0, T = 1) channel, 
• • ■ a bound state, the deuteron, is formed. Two-nucleon states in nuclear matter are strongly modified when increasing 
the density. In particular, the bound state will disappear at the so-called Mott density ||||. 

Another aspect of the nucleon-nucleon interaction is the formation of quantum condensates. At low temperatures it 
is well known that in nuclei, nuclear matter, and neutron matter (neutron stars) superfluidity can arise in the isospin 
triplet channel HJ^o' Nucleonic pairing is a well-known effect in the structure of nuclei. In the isospin triplet channel 
the influence of proton-proton (pp) or neutron-neutron (nri) pairing on the binding energy of nuclei has extensively 
been investigated, cf. Q. 

However, under certain conditions the interaction in the proton-neutron (pn) isospin singlet channel may be even 
stronger. For instance, in the two-particle system a bound state, the deuteron, arises. In nuclear matter, at densities 
below the Mott density a bound state (quasi-deuteron) can exist so that at low temperature the Bose-Einstein 
condensation of these quasi-deuterons may occur. An interesting feature of the isospin singlet pairing in symmetric 
nuclear matter is the cross-over from Bose-Einstein condensation of deuterons at low densities to BCS neutron-proton 
pairing at high densities p|,|lO|]. 

In spite of the relatively strong interaction, triplet pairing seems less apparent in nuclear structure systematics (see 
the studies of Goodman on T = pairing in nuclei pj). However, it should be essential for nuclei near the symmetry 
line N » Z. This phenomenon becomes of importance for heavier N as Z nuclei p2| as will be produced in the new 
radioactive beam facilities. 

A further interesting effect is the possible occurrence of higher order condensates such as a- like quartetting |13L[l4j . 
At present, the signatures of isospin singlet pairing and the relation to quartetting have not so clearly been worked out. 
An interesting point is to identify signatures of quartetting in finite nuclei, see UM . One of the possible consequences 



is the contribution of isospin singlet pairing and four-particle correlations to the binding energy of nuclei, that will be 
discussed below. 

The treatment of correlations in the many-nucleon system is more simple for infinite nuclear matter. The new 
results with respect to isospin singlet pairing and quartetting |15| should be of interest for the evaluation of binding 
energies of finite nuclei as well. One possible way to apply results of infinite matter calculations to finite systems is 
the local density approximation (LDA) which was elaborated successfully in atomic, molecular and condensed matter 
systems, sec [ |16| and references therein quoted. We will use this approach to perform exploratory calculations in 
finite nuclei. Compared with shell-model calculations, the LDA is by far more simple and has been proved to give 
adequate results in the case of pp and nn pairing |17j . We will use this approach to give estimates for the effects of 
isospin singlet pairing and quartetting on nuclear binding energies. 

The question of different kinds of pairing in proton-rich nuclei has been investigated by various authors. In addition 
to the T = \ pairing, T = np-pairing modes have been investigated to explain the Wigner energy in N = Z nuclei 
p8| . A more extended discussion of the origin of the Wigner term is given in |l9| where a method has been developed 
to extract the Wigner term from experimental data. Both, empirical arguments and shell-model calculations suggest 
that the Wigner term can be traced back to the isospin T = part of nuclear interaction. However, it has also 
been found that the Wigner term cannot be solely explained in terms of correlations between the neutron-proton 
J = 1, T = (deuteron-like) pairs. One should, however, realize that in a finite nucleus the intrinsic quantum 
numbers of the deuteron can be mixed with the orbital motion of the nucleons in the shell model potential. Recently, 
pairing and the structure of the pj-shell N rs Z nuclei have been discussed in J20| . There, the isovector and isoscalar 
pairing interaction has been studied. It has been found that the Wigner energy cannot purely be explained as a 
pairing effect, considering only the zero angular momentum (L = 0) channel. 

Calculations of the pairing and isospin symmetry in proton-rich nuclei were performed [EM. Near N w Z, a steep 
decrease of the isoscalar proton- neutron pairing energy is found with increasing \N — Z\. An interesting fact is that in 
50 Cr proton- proton and neutron- neutron pairing is more reduced with temperature than the isoscalar proton-neutron 
pairing ||. 

Proton- neutron pairing has been discussed also with respect to other properties, see |lq| . For instance, a cranked 
shell model for the description of rotational bands in N w Z nuclei has been formulated in |23[ taking the isovector 
proton- neutron pairing explicitly into account. More detailed discussions with respect to the binding energies can be 
found in fll4J . Proton-neutron pairing and its relation to isospin symmetry has been studied within the BCS theory 
in [p4| . Recently, the role of the proton-neutron interaction in N m Z nuclei and its consequences for pn-pairing has 
been investigated in p5[ . 

We will proceed as follows: In Section 2 we analyze the binding energies of nuclei near N — Z. The basic formalism 
of LDA is given in Section 3. We proceed then with the presentation of numerical results in Section 4 and draw some 
conclusions in Section 5. 

II. BINDING ENERGIES OF NUCLEI NEAR N = Z 
A. Filters for isospin singlet condensates 

The empirical binding energies of nuclei are well approximated by the Bethe-Weizsacker formula, containing the 
contributions of bulk, surface, Coulomb and asymmetry energy. Additional contributions originate from pairing in 
the isospin triplet channel (pp, nn) and shell structure effects, for instance the behavior near magic numbers. 

We assume that the binding energy B(Z. N) can be decomposed into different contributions as 

B(Z, N) = B hulk (Z, N) + B suli (Z, N) + B Coul {Z, N) + B, symm (Z, N) 

+ B shcll (Z, N) + B pair (Z, N) + AB(Z, N) . (1) 

The term AB(Z,N) contains the effects of correlations between protons and neutrons (isospin singlet (pn) pairing, 
quartetting) which are of interest here. In nuclear matter, the occurrence of different condensates has been intensively 
investigated. Isospin singlet (pn) pairing is strong in symmetric nuclear matter for densities up to the saturation 
density. The dependence of the pn pairing gap on the asymmetry of nuclear matter for different densities has been 
investigated in [£6|,g7|]. If the densities of protons and neutrons, or the corresponding chemical potentials fj, p , fj, n are 
sufficiently different, a condensate cannot be formed. 

Analogously, the effects of isospin singlet (pn) pairing and quartetting are expected to be strong near N = Z 1 
and disappear if the absolute value of the difference N — Z becomes large. Therefore, it seems to be reasonable to 
investigate the behavior of 



AB(Z,N) = b N _ z (Z) (2) 

in dependence of TV — Z, considering Z as a parameter. 

To select out the contribution AB(Z, TV) to the binding energy, filters can be applied which eliminate the other 
contributions given in (|l|) to a large extent. We consider the following (horizontal and vertical) filters 

h{Z, TV) = 2B{Z, N) - B{Z, N - 2) - B(Z, TV + 2) 

-2B(Z - 2, TV) + B(Z - 2, TV - 2) + B(Z - 2, TV + 2) , (3) 

v(Z, TV) = 2B(Z, TV) - B(Z - 2, TV) - B(Z + 2, TV) 

-2B(Z, N + 2) + B(Z - 2, TV + 2) + B(Z + 2, TV + 2) . (4) 

These simple filters contain the difference of second order differences. Therefore, smooth dependences on Z, TV nearly 
cancel out, up to second order. This concerns not only the contributions of bulk, surface, and Coulomb energy, 
but also the asymmetry term. The isospin triplet pairing effects compensate because only nuclei of equal parity are 
considered. Similarly, shell effects compensate to a certain extent assuming that they are related only to the respective 
values of the numbers Z or TV. 

From these primitive filters (ph, (0), other filters can be constructed by superposition. The filter used in [H to 
evaluate the Wigner energy can be constructed as 

™~;-&4-»H»(4.4-»). 



or, more generally 



Another, more symmetric filter 



w{Z,N) = -\v{Z,N-2)-\h{Z,N-2). (6) 



g(Z,N) = ±h(Z,N)-±h(Z + 2,N) (7) 

was considered in |2S] ]. 

We give the results for the filter h(Z, TV) in Fig. |l|. The dependence of h(Z, TV) on Z is shown for different parameter 
values TV — Z, where even and odd values of Z are shown separately. It is clearly seen that the filter h(Z, TV) becomes 
small for large TV — Z . Furthermore, the average value decreases with increasing Z . Also, even-even and odd-odd 
numbers of Z, TV show a different behavior. Similar results can be obtained using the filters (0), (ph, or (R). 

The relation to the quantities bi(Z) introduced above can immediately be given. For the filter h(Z, TV) we find 

h(Z, TV) = -b N -z-2(Z) + 2b N _ z (Z) + b N - Z {Z - 2) 

-b N _ z+2 {Z) - 2b N _ z+2 (Z - 2) + b N _ z+i {Z - 2) . (8) 

Since only differences of the quantities bi(Z) are determined by the filter h(Z,N), their values can be reconstructed 
using additional assumptions. To simplify the extraction of the quantities bi(Z) we assume that bz-N^Z) is strongly 
depending on the asymmetry (Z — TV), but only weakly depending on the mass number (Z), after eliminating a possible 
even-odd staggering. Therefore, we approximate bi(Z — 2) rs bi(Z) so that h(Z,N) ta —bN-z-2(Z) + 3bN-z{Z) — 
3b]y_z+2(Z) + bfli_z+4:{Z)- According to the increase of asymmetry discussed above, we assume that the bi(Z) 
disappear for large absolute values of i. Taking bi = for |T| > 6, the results for the quantities bi(Z) are shown in 
Fig-! 

A first important result is that the quantities bi(Z) depend mainly on the absolute value of i. As a consequence, 
h(Z, Z — 2) nearly coincides with —h(Z, Z), whereas h(Z, Z — 1) is close to zero, in contrast to h(Z, Z + 1). In detail, 
the values for bo, b 2 , 64 were found from the solution of equations containing h(Z, Z), h(Z, Z + 2), h(Z, Z + 4), whereas 
the values for b\, 63, 65 were found from the solution of equations containing h(Z, Z + 1), h(Z, Z + 3), h(Z, Z + 5). 
To derive the average behavior with respect to Z, we considered the ratios bi/b as shown in Fig. [| The average of 
these values over the whole range of experimentally accessible data (8 < Z < 30) are given in Fig. [|. They show a 
decrease with increasing |i| = |TV — Z\, so that 65 is nearly zero. 



B. Further filters 

A further interesting property of the quantities bi(Z) is their dependence on Z. First we are interested in a possible 
even-odd staggering. This can be deduced from the values shown in Fig. g. Alternatively, we also can consider, as a 
special indicator to extract even-odd staggering, the filter 

c(Z, N) = (-if \B(Z, N) - B(Z, N + 2)- B{Z + 2, N) + B{Z + 2,N + 2) 

+ +B(Z + 1, N - 1) - B(Z + 1, N + 1) - B(Z - 1, N - 1) + B(Z - 1, N + 1)] . (9) 

We decompose 

bi(Z) = HZ) + (-ifSb^Z) (10) 

and assume, that the remaining dependence on Z is smooth, as already used above. Then we have 

c(Z, N) = -2Sb N - Z -2(Z) + 48b N -z(Z) - 25b N - Z+2 (Z) . (11) 

Now, we can perform the reconstruction of Sbi(Z) along the lines given for bi(Z). The filter c(Z, N) is shown in Fig. 0, 
the parameters Sbi(z) in Fig. |as a function of Z, and the averaged values are shown in Fig. |7] vs the asymmetry 
i — N — Z. We see that also there the staggering contribution disappears with increasing \i\, however, the statistical 
errors are large. The global dependence of bi(Z), 5bi(Z) on Z cannot easily be extracted. After a steep decrease for 
small Z a flattening is observed at higher values of Z, see also Figs. |[ |[ 

From other approaches, see ]19fl , the Wigner energy Ew is introduced which occurs in the additional binding due 
to the np-pair correlations 

B n p,paii = ~ e np\A)Tl n p + Ew (1^) 

with 7T„ p = (1 — (— 1) JV )(1 — (— l) z )/4. The first contribution to the np pairing energy (e np ) represents additional 
binding due to the residual interaction between the two odd nucleons in an odd-odd nucleus. The Wigner energy Ew 
is believed to represent the energy of collective np-pairing correlations. It can be decomposed into two parts: 

E w = W(A)\N-Z\ + d(A)Tr np 6 NZ (13) 

Filter used in Q can be constructed as given by (||) and 

,/.n 1 (A A \ 1 (A A\ 
d ^ = -2 C {2- 3 ^- 1 ) + 2 C {-2- 2 '2)- (W) 

The d-teim represents a correction for N — Z odd-odd nuclei. Estimates suggest that the ratio d/W is constant, values 
1 and 0.56 have been reported |lq| . Both filters can be expressed by superposition in terms of the more simple filters 
h(Z,N), v(Z,N), c(Z, N) given here. A more general filter is introduced as d{Z,N) = c(Z- 3, N-l)/2 + c(Z -2, Z) 
as an average over neighboring values of c. Some values are shown in Fig. H. 

The results of our phenomenological treatment can be summarized as follows (some preliminary results have been 
presented in p9||): 

(i) There is a contribution bi(Z) to the nucleon binding energy which is due to proton- neutron correlations of the 
type of a isospin-singlet pairing or a quartetting condensate. 

(ii) This contribution depends on the absolute value of the asymmetry parameter \N — Z\. It has a maximum 
magnitude for symmetric nuclei N = Z and decreases with increasing asymmetry, disappearing near \N— Z\ = 4. 

(iii) It shows an even-odd staggering as function of Z. 

(iv) On the average, it decreases with increasing Z, steep for small values of Z, but fiat for large Z. 



III. LOCAL DENSITY APPROXIMATION 

A theoretical interpretation of the contribution to the binding energy due to pn pairing and quartetting could be 
given by the Local Density Approximation (LDA). In contrast to shell-model calculations, the LDA is by far more 
simple and has been proved to give adequate results in the case of pp and nn pairing [[17J . 

As well known from quantum statistics of the inhomogeneous fermion gas, the energy and wave function of the 
ground state of a many-fermion system can be calculated within a variational approach. The energy density is 
considered as a functional of the fermion density, which, in the case of a nucleonic system, depends in addition to 
space coordinates also on spin and isospin. It can be decomposed into kinetic, potential, exchange and correlation 
energy. Within a gradient expansion, in lowest order the energy density depends only on the local values of the 
nucleon density. As a consequence, the exchange and correlation energy can be approximated using results from 
nuclear matter calculations. 

In particular, the contribution to the exchange and correlation part of the energy due to the formation of a 
condensate can be evaluated within nuclear matter theory. There is an extended literature on isospin triplet (pp, nn) 
pairing. More recently, also isospin singlet (pn) pairing has been considered 0], which for symmetric nuclear matter 
may become stronger compared with isospin triplet pairing at subnuclear densities because of the more attractive 
nucleon-nucleon interaction in the isospin singlet channel. 

A. Pairing vs. quartetting in symmetric matter 

A standard way to describe quantum condensates in many-body systems is the method of thermodynamic Green 
functions. Treating the two-particle Green function in ladder Hartree-Fock approximation, an effective wave equation 
(in matrix notation) ip\ — K2(E\) ip\ for the quantum state A can be derived. Explicitly this reads 

^(12) = ^A- 2 (12,l'2',e A )^ A (l'2') (15) 

with 



1'2' 



K 2 (12, 1'2', z) = 1/(12, 1'g) \ J^ _ ^ . (16) 



The influence of the medium is contained in the single-particle energy 

E 



e(l)=p?/2m + yV(12,12)/(2) (17) 



and in the Pauli blocking term [1 - /(l) - /(2)]. Here, /(l) = [exp{e(l)/T - fi/T} + l]" 1 is the Fermi distribution 
function and '1' denotes momentum, spin, and isospin coordinates, whereas T^(12, 1'2') is the antisymmetrized matrix 
element of the two-body interaction. 

The transition to a superfluid state is obtained from the Thouless criterion as described by the Gorkov equation 
tp2 = A^(mi + M2) ip2- Depending on the respective channels considered, it allows the determination of the critical 
temperatures T° or T t c for the isospin singlet and triplet channels, respectively, as a function of the chemical potential. 

The solution of the Gorkov equation has been considered by different authors using realistic bare nucleon-nucleon 
interactions. It has been found that in comparison with the isospin triplet channel, in the isospin singlet channel 
the transition to superfluidity should arise at relatively high temperatures [P Pqj27[| , see also Figs. |], |lO]. This is a 
consequence of the stronger interaction in the isospin singlet channel which leads to the formation of the deuteron in 
the low-density limit where / <C 1. Estimates give a value of the critical temperature up to T s c ss 5 MeV at one third 
of the nuclear matter density. At the same time, at zero temperature a large gap arises M. 

In a recent letter |15| it has been shown that in a certain region of density, pairing has to compete with quartetting. 
It has been found that under certain conditions in symmetric nuclear matter the transition to isospin singlet pairing, 
which is stronger than triplet pairing, will not occur because the quartetting transition occurs before that. Within 
a cluster-mean field approach [ppQ], the critical temperature for the quartetting transition was obtained from the 
equation 

G 4 (1234 > 1'2'3'4' ) ,)= WWWW> *^-*«" 



34(1234) z - e 4 (1234) 






1"2"3"4" 



#4(1234, 1"2"3"4", z)G 4 (l"2"3"4", l'2'3'4', z) , (18) 



#4(1234, 1'2'3'4', z) = V(12, 1'2') 



/(l)/(2) <WW 

0j(12) z-e 4 (1234) 



perm. . 



(19) 



where we use the abbreviation e„(12 . . . n) = e(l)+e(2) + - ■ - + e(n), and g n (12 . . ,n) = [exp(e„(12 . . .n) — rifj,)/T— 1] _1 
being the Bose distribution function. The instantaneous part of interaction kernel is obtained by using the technique 
of Matsubara Green functions as where the terms obtained by renumbering are not given explicitly. We have used 
the identity /(l)/(2) • ■ • f(n) - /(l)/(2) ; ■ • /(») = g-\l2 . . . n)/(l)/(2) • • • f(n) with / = 1 - /. The solution of 
the equation ip^ = K^Afijip^ gives the critical temperature for the onset of quartetting as a function of the chemical 
potential as shown in Fig. 0, or the density as shown in Fig. [lO. Within an estimate by using a variational calculation, 
the transition to quartetting beats the transition to isospin singlet pairing if the density is smaller than 0.03 fm _ , 
see Fig. [io. 



B. Gap equation and condensation energy for asymmetric nuclear matter 

For infinite nuclear matter, the gap energy at zero temperature as well as at finite temperature has been investigated 
for pairing in the different channels in dependence on nucleon density and isospin asymmetry. In particular, it has 
been found that the gap energy in the isospin singlet channel is strongly reduced for increasing asymmetry, and the 
transition to superfluidity is possible only for asymmetry values a — (n n — n p )/(n n + n p ) < 0.35 [E6 27|. Also, the 
critical temperature is strongly suppressed with increasing asymmetry as it can be calculated from the gap energy as 
well as directly from the solution of the BCS equation. 

We give some relevant expressions for a fermion system interacting via a separable potential 



V rr , = -A J] w{k)w{k') a\{P/2 + k) a\,{Pj2 - k) a T ,(P/2 - k') a T (P/2 + k') 



(20) 



P,k,k' 



We use a Yamaguchi type of potential |3l| with w{k) = (k 2 + k 2 ) -1 , and n — 1.4488 fm . The interaction strength 
in the 1 So channel is only about 70 percent of the strength in the 3 Si channel which in its original form pl| is chosen 
to reproduce the deuteron binding energy as well as the low-energy behavior of the free scattering phase shifts. We 
perform an exploratory calculation and consider the interaction strength as a parameter that will be adjusted below. 
Furthermore, only zero angular momentum is considered. Separable representations of more realistic interactions can 
be found in the literature p2] |. 

The interaction is treated in mean-field (Hartree-Fock-Bogoliubov) approximation, allowing for an isospin-singlct 
pair amplitude at zero total momentum P. Diagonalizing H MF — [i p N p — /z n iV n using the Bogoliubov transformation, 
we obtain the gap equation 



A(fc) = Xw{k)J2^( k ') 



A(fc') 



V / fe(fc') + Cn(fc')) 2 +4A2(fc') 



[l-f(E+)-f(E-] 



with 



E 



v (&(*) + Uk)) 2 + 4A2(fc) ± fa(k) - Uk)) 



(21) 



(22) 



and £, p (k) — e p (k) — /i p , £ n (fc) = e n (k) — Un, f(E) — (e E ' T + 1) — 1 . t T (k) are the single-particle energies including the 
shift due to a mean field as given above dl^). From the self-consistent solution of the gap equation, besides the trivial 
solution A(fc) = also a solution A(fc) = g w(k) with a finite value of g may occur. 

The shift in the energy density due to the formation of a gap (condensation energy density) follows as J33| 



L\C\Tl p: Tl n J — |£pair \TT> pj Tl n J c-norm \^p-j ^n )\ 



-y:\- 

V L^i 1 o 



V 



1- 



Uk)+Uk) 



Vfe( fc )+^(fc)) 2 +4A2(fc) 
A 2 (fc) 



(e p (k) + e n {k)) 



y/(Uk)+Uk)) 2 +4A 2 (k) 



[1 + f{Et) + f(E k 



-^(K) 



{e P (k) - e„(fc)) + (e p (fc) + e n (k)) 



U k )+Uk) 



Vfe(fc)+CnW) 2 +4A 2 (fc) 



-J(Ek) 



A 2 (fc) 



(e n (k) - e p (k)) + (e p (k) + e„(fc)) 



A 2 (fc) 



&(*)+&(*) 



Vfe(fc)+e«(fc)) 2 + 4A2(fc) 



-£ P (k)f(e p (k) - n p ) - e n (k)f(e n (k) - ^ n ) 



(23) 



The chemical potentials are given by the normalization to the densities of the corresponding nuclcons 

2 -r-^ 

yL/( £ rW -Mr) =n r . 



(24) 



As usual, it is assumed that the normalization condition in the paired state gives no essential change in the corre- 
sponding chemical potentials Q . 



C. Finite nuclei density profiles 

For infinite nuclear matter, the energy density is calculated for homogeneous densities n T . In finite nuclei, the 
densities n T {r) are depending on position r. Then also the quantities considered above which are functions of the 
densities now are parametrically depending on r. 

If the density distribution of protons and neutrons is known, the gain of the binding energy due to np pairing 
(condensation energy) can be estimated in LDA by the integral 



B np « Ait 



r A£(n p , n„; r)dr 



(25) 



For nuclei with N neutrons and Z protons and A = N + Z, we have to determine the density profiles of protons 
and neutrons. Exploratory calculations will be performed taking the nucleon densities from a simple potential model, 
normalized to the corresponding numbers of protons or neutrons. 

The nucleons feel a mean field phenomenologically defined as 



V p (r) = V^\r) + V cou \r), V n (r) = ^V) • 

We adopt the Shlomo parameterization for the mean field fl34|] : 

V? 



rynucl 



1 + exp((r - R)/d) ' 



V° = -V°+ T V sym (N -Z)/A, 



(26) 



(27) 



with the parameter values V° = 54 MeV, V sym = 33 MeV, and d = 0.7 fm. The radius of the nuclear potential is 
given by the implicit equation 



R = 



1.12A 1 / 3 + 1.0 

[l + (^/i?) 2 ]V3 

For the Coulomb potential we use the charged sphere formula 



T/XOul 



Ze 2 
2R C 



3- — 
6 Rl 



Ze 1 



Q(R C -r) + 0(r - R c ) , 



where e 2 = 1.44 MeV fm, Q(x) denotes the step function. The Coulomb radius is given by 

2 _ 5 ,_ ^ 2 l + 10/3(^/C) 2 + 7/3(Wg) 4 
Kc - 3 < r >- L 1 + (Trz/Cr ' 



(28) 



(29) 



(30) 



c 1 - 12A1/3 m) 

[1 + (ttz/C) 2 } 1 / 3 ' l ' 

where < r 2 > is m.s. radius from the nuclear charge density, z — 0.54 fm. 

Within the Thomas-Fermi approximation, at zero temperature the local density is given by 



Mr) = ^(fcf) 3 6(A r -W(r)) 



1 

2m . , . . 

— (X T -V T (r)) 



1/2 

(32) 



K 

where V T (r) is the total potential (|2q). The chemical potentials A T are determined by the constraints 

d 3 r n p {r) , N = f d 3 r n n (r) . (33) 

Obviously this simple treatment gives only a first estimate of the nucleonic densities. They are correctly normalized 
and account also for the Coulomb repulsion, but can be improved by considering the quasiparticle dispersion relation 
or shell effects in the nuclear density. 

IV. NUMERICAL RESULTS 
A. Condensation energy 

In the LDA, we proceed in the following way. First, we calculate the density profiles n T (r) for protons and neutrons 



(see Sec. [II C). As typical examples the density profiles for °Ar and 40 Ti are shown in Fig. 11. Then, with the 
total local density and the asymmetry [n n (r) — n p (r)]/[n n (r) + n p (r)\ as inputs the local gap function A(fc;r) is 
obtained from the solution of the gap equation (Ell) . This is performed within a self-consistent Hartree-Fock scheme, 
which gives the local single-particle energies e T (fc;_r) depending on the isospin variable. Having the pairing gap at our 
disposal the local condensation energy density p3) is determined. Finally, integration over the whole nucleus ( J25J ) 
gives the contribution to the binding energy due to np pairing. 

The coupling strength A of the separable interaction is considered as a parameter which should be taken as a 
phenomenological quantity. To reproduce the average Wigner energy in the mass number region 20 < A < 100 
(2.62 MeV), we have taken the value A = 92.35 MeV fm 3 in the % channel and A = 131.50 MeV fm 3 in the 3 Si 
channel. Usually the deuteron binding energy is used to adjust this parameter. However, within our simple model we 
should take into consideration that the interaction (|20|) is an effective description, simplifying different contributions 
as the different channels, short distance repulsion, spin-orbit coupling etc. Furthermore, some suppression of T = 
pairing in symmetric nuclear matter due to medium polarization is expected. The fact that the np pair is bound 
(deuteron) whereas nn and pp pairs are not is essentially due to the tensor force leading to the <i-wave component in 
the deuteron. Without this component the pn interaction in the T = channel would hardly be different from the nn 
or pp interactions. The fate of the tensor force in the nuclear medium is, however, a much debated subject in nuclear 
physics and it is quite possible that the tensor force is much more screened than the other parts of the nuclear force 
p5[ . In this sense the use of a bare interaction in the pn (T — 0) channel may be more questionable than it is in the 
T = 1 channel. 

The isospin singlet gap for symmetric nuclear matter is shown in Fig. |l2j for different temperatures. Below normal 
nuclear matter density, which is of relevance here, the difference between T = and T — 0.5 MeV temperature is 
small. The dependence of the isospin singlet gap on the nuclear matter density is shown in Fig. [13j for the temperature 
T = 0.75 MeV and different asymmetries. Thegap is strongly reduced for increasing asymmetries and temperatures 
confirming the results obtained previously in |26j,|27J . 

For a nucleus with A = 40, the average gap on the Fermi level |17J] 



A = Y^ [ d 3 rA(k F (ry,r)[n T (r)}^ /j2 f d ^ 



Mr)] a (34) 



is shown as function of the asymmetry in Fig. [tl] for different temperatures. 

The condensation energy as a function of asymmetry is shown in Fig. [fq for different temperatures. For the 
calculation, a fixed A — 40 was assumed. Below T = 0.5 MeV, the dependence on temperature is negligible (see 



also Fig. |lj), but becomes strong for T > 1 MeV. The influence of the Coulomb interaction taken into account for 
the calculation of the density profiles destroys the symmetry with respect to N — Z = 0. The Coulomb effect is to 
increase the overlap between neutron and proton densities in the tail of the density profiles, as show in Fig. [[T| As 
a consequence the pairing gap is slightly enhanced. Furthermore, the steep decrease of the condensation energy near 
N — Z = 5 is also shown. This is in correspondence to the findings given in Fig. 0. 

B. Quartetting 

The additional contribution due to quartetting seems to be high in the region of light nuclei, where the a cluster 
model is a good approximation. The strong even-odd staggering is reduced at higher masses. 

To give an estimation of the effect of quartetting, the evaluation of the np condensation energy ( |25| ) was repeated 
for symmetric nuclei (Z — N), where the magnitude of the gap A(fc;r) was increased in the density region where 
quartetting can occur, i.e. at densities below 0.03 fm _ . In detail, the gap was increased by a factor, which was 
obtained from the ratio T^/T^ of the corresponding critical temperatures at given density n, as shown in Fig. |l0| . 

Comparing with calculations neglecting quartetting, the gain of binding energy due to quartetting has been evaluated 
for different nuclei with Z = N. Exploratory calculations for nuclei of medium size (A w 100) show that the 
contribution due to quartetting is almost zero but may become large for small A. For instance, the calculation for 
12 C and 16 give an additional contribution to the condensation energy due to quartetting of 10.3 and 9.6 percent, 
respectively, if compared with isospin singlet pairing. 

C. Excited Nuclei 

After discussing the contribution of pn-pairing to the nucleon binding energy, comparing with the pairing energy 
in asymmetric nuclear matter, it is of interest to discuss also the effect of excitations. In nuclear matter, excitations 
are well understood in the context of finite temperatures. It is expected that the effects of condensates arc decreasing 
with increasing excitation. 

To investigate the effect of excitations on the formation of condensates in finite nuclei, we analyse the 2 + excitations 
of even-even nuclei. In particular, we used the w-filter for the analysis and compared u>2+(Z, Z) for the excited nuclei 
with wq(Z, Z) for the ground state nuclei, see Fig. 061 We determined the mean value of the ratio from Z = 10 to 
Z = 26 and obtained w 2 + (Z, Z)/w {Z, Z) = 0.64 ± 0.19. 

This result can be compared with the influence of finite temperature on the pairing in nuclear matter. The average 
2 + excitation energy of even-even nuclei, taken for the interval 10 < Z < 30, is 1750 keV. Performing a finite 
temperature Thomas-Fermi calculation for a nucleus of medium proton number Z — 20, this excitation energy would 
correspond to a temperature of about 1 MeV. 

As shown in Fig. [jjj, similar to the decrease of the gap we observe also a decrease of the value bo(Z) with increasing 
temperature. The result of the calculation is in agreement with the empirical value given above. 

V. CONCLUSIONS 

In this work we have shown that isospin-singlet pairing and a like quartetting may contribute to the binding energy 
of nuclei very close to the symmetry line N = Z. These contributions are relatively large for smaller nuclei (Z < 20). 
For medium mass nuclei neutron-proton pairing in the isoscalar channel disappears already for \N — Z\ = 4. This 
stems from the very rapid decrease of the isoscalar pairing as a function of the unbalance in the Fermi energies of 
protons and neutrons. These facts can explain the origin of the Wigner term in the mass formula as well as the 
empirically determined AB(Z,N) or bz-N(Z), respectively (see Eqs. (Q), (|l0|)). An enhancement of the isoscalar 
pairing contribution to the binding energy is obtained if in addition quartetting is taken into account. An interesting 
effect is that the reduction of the condensation energy with increasing excitation of the nuclei seems to be in agreement 
with empirical data. 

Our calculations are exploratory in the sense that they were performed in the rather crude LDA approach. However, 
please notice that the LDA has yielded in the past quite reasonable results on the average, i.e. for a gap averaged 
over the shell effects JlT]]. Therefore, we think that our results give a quite reliable first orientation of the effect. 
The approach should be improved in several respects. First a more realistic force should be employed. Second shell 
effects must properly be included. Eventually, effects of number projection, even-odd staggering and pair fluctuations 
should also be investigated. Such studies shall be performed in the future. It is the hope that isoscalar pairing and 



quartetting will give us precious hints on the effective neutron-proton interaction in a nuclear medium as well as very 
interesting clustering and condensation phenomena in nuclei. 
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FIG. 1. The filter h(Z, Z + i), Eq. (p|), is given for various i as it can be extracted from the experimental binding energies. 
Results are separately shown for even-even (a), odd-odd (b), and even-odd/odd-even (c) nuclei. 
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FIG. 2. The parameter bN-z(Z) for even and odd proton number Z as it is derived from the filter h. 




FIG. 3. A selection of several ratios bi/bj as a function of Z. The dashed lines represent the average values of the ratios. 
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FIG. 4. Averages of the ratio 6jv_z/bo as a function of the absolute value of the difference between protons and neutrons. 
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FIG. 5. The filter c(Z, Z -\-i) for various i. Results are separately shown for even-even/odd-odd (a) and even-odd/odd-even 
(b) nuclei. 
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FIG. 6. The parameter SbN-z(Z) as derived from the filter c, Eq. fflf), for even and odd proton number Z. 
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FIG. 7. Averages of the ratio SbN-z/Sbo as a function of the absolute value of the proton-neutron difference. 
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FIG. 8. For demonstration purposes the filter d(Z, Z + i) for various i is given considering only even-even and odd-odd nuclei. 
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mation for 40 Ar and 40 Ti at zero temperature. 
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FIG. 16. For the filter w(Z, Z) the ratio w 2 +/ioo of 2 + excited nuclei over ground state nuclei is given for even proton and 
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